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Q ! We report on two Chandra observations of the 3-Myr-old pulsar B1929+10, 
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which reveal a faint compact (~ 9" x 5") nebula elongated in the direction per- 
pendicular to the pulsar's proper motion, two patchy wings, and a possible short 
(~ 3") jet emerging from the pulsar. In addition, we detect a tail extending up 

^ \ to at least 4' in the direction opposite to the pulsar's proper motion, aligned 

with the ~ 15'-long tail detected in ROSAT and XMM-Newton observations. 
The overall morphology of the nebula suggests that the shocked pulsar wind is 

'Ec ! confined by the ram pressure due to the pulsar's supersonic speed. The shape of 

the compact nebula in the immediate vicinity of the pulsar seems to be consistent 

^ I with the current MHD models. However, since these models do not account yet 

K*" I for the change of the flow velocity at larger distances from the pulsar, they are 

[~^ I not able to constrain the extent of the long pulsar tail. The luminosity of the 

whole nebula as seen by Chandra is LpwN ~ 10^° ergs s~^ in the 0.3-8 keV band, 
for the distance of 361 pc. Using the Chandra and XMM-Newton data, we found 
that the pulsar spectrum is comprised of non-thermal (magnetospheric) and ther- 

Q I mal components. The non-thermal component can be described by a power-law 

model with photon index F ^ 1.7 and luminosity -Z^psr*^ ~ 1-7 x 10^° ergs s~^ 

k>( ; in the 0.3-10 keV band. The blackbody fit for the thermal component, which 

H I presumably emerges from hot polar caps, gives the temperature kT ?» 0.3 keV 

and projected emitting area A± ~ 3 x 10^ m^, corresponding to the bolometric 
luminosity Lboi ~ (1^2) x 10^° ergs s"^. 

Subject headings: pulsars: individual (PSR B1929-I-10) — stars: neutron — X- 
rays: stars 



1. Introduction 

Pulsar wind nebulae (PWNe), generated when the magnetized relati vistic pulsar winds 



shock in the ambient medium, have been observed around ~ 50 pulsars (see lKargaltsev fc Pavlov 
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20081 . for a recent review). The shocked pulsar winds produce synchrotron radiation from 
radio frequencies through 7-rays, reveahng properties of the winds and the ambient medium. 

The PWN morphology depends on the pulsar's velocity, pressure, temperature, and 
magnetic field of the ambient medium, and the geometry of the wind outflow. For instance, 
the torus-jet PWN structures obser ved around young pul s ars in supernova re mnants (SNRs), 
such as the Crab and Vela PWNe ( jWeisskopf et al.ll2000l : iPavlov et al.ll2003l ). are commonly 
interpreted as anisotropic outflows with equatorial and polar components confined by the 
pressure of the hot gas in the host SNR interiors, where the speed of sound exceeds the 



pulsar's speed ( iKomissarov &: LyubarskyI l2003l : iDel Zanna et al 



2OO4J ). When the pulsar is 



moving wit h a supersonic speed with respect to the amb ient medium, bow-shock PWNe 



are formed (JBucciantini et al.ll2005l: iRomanova et al 



J174 7-2958 (the Mouse PWN; lOaensler et al.l l2004h and PSR B1957-20 (IStappers et al. 



20051 ). such as the PW Ne around PSR 



2OO3I ). The pulsar's supersonic speed causes the ram pressure to exce ed the ambient pressure , 
so that the wind termination shock acquires a bullet-like shape (e.g. iBucciantini et al.ll2005l ) 
and forms a PWN, sometimes with a long tail extending behind the pulsar. The pulsar's 
speed becomes supersonic as it encounters a relatively cold interstellar medium (ISM), either 
after leaving the SNR in which it was born, or when the host SNR dissolves in the ISM 
background. These last stages of the SNR evolution occur after ~ 10^ years. 

The majority of middle-aged (~ 10^-10^ yr) and old (> 10^ yr) pulsars are expected to 
move supersonically through the surrounding ISM, and as long as their winds are sufficiently 
strong, they are expected to produce observable bow-shock PWNe. Observations of these 
objects provide an opportunity to study bow-shock morphologies and probe pulsar winds 
and their surroundings at later stages of pulsar evolution. 

Until now, only a handful of PWNe ar ound middle-aged pulsars have been det ected: 
e.g., PSR J0538+2817 [ ' 

Geminga (|Caraveo et al. 



quencies (JGaensler et al 



2004 



Romani fc Nglboosh , PSR B0355+54 JMcGowan et al.ll2007h . and 



Pavlov et all l2006f ) in X-rays, PSRB09 06-49 at radio fre - 



19981 



I, and PSRB1951+32 in X-rays and Ha JMoou et al.ll2004h . 



However, one of the first bow-shock PWN candidates was detected serendipitously around 
a much older p ulsar, PSR B1929+10 ( hereafter B1929). This is one of the nearest pulsars 
{d = 361^8° pc; IChatterjee et al.ll2004l ). with a period P = 226.5 ms and period derivative 
P = 1.16 X 10^^^ s s~^. Despite the large spindown age, r = P/2P = 3.1 Myr, its spin- 



down power, E = Att'^IP/P^ = 3.9 x 10^'^ ergs s ^ for the moment of inertia / = 1 x 10^^ 



is still relativ ely high. As the pulsa r's transverse velocity in the plane of the sky, 
177^5 km s~^ ( IChatterjee et al.l 120041 ). substantially exceeds the typical ISM sound 



g cm 

V^ = 

speed, one can expect that the puls ar wind outflow for ms a bow-shock PWN. An indication 

of such a PWN was flrst noticed by lWang et al.l (119931 ). who detected a very long (~ 1.6 pc) 
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tail-like structure behind this pulsar in ROSAT data. The tail behind B1929 had been the 
longest extended structure associated with a compact Galactic object until the recent dis- 
co yery of a G-pc-lon g X-ra y tail behind the middle-ag ed pulsar PSR J1509— 58 50, reported 
by iKargaltsev et al.l (120061 ) and described in detail by iKargaltsev et al.l (1200811 Similar (al - 
beit shorter) tails have been found behind a few other pulsars ( iKargaltsev fc Pavlovl l2008l ) 
suggesting that such extended PWN morphologies might be ubiquitous. 

The B1929 pulsar has been extensively observed in X-rays, in an attempt to study the 
magnetospheric and thermal components of its emission and compare them with the models 
for pulsar radiation. In particular, the thermal emission from polar caps heated by relativistic 
parti cles accelerated in the pulsar mag netosphere is predicted by the current pulsar models 
(e.g. jHarding fc Muslimovll200ll . l2002l ). but a limited sample of the observed old pulsars sug- 



gests that the nonthermal magnetospheric emission with a power-law spectrum dominates 
at higher X-ray ene rgies, E >,2 keV, while the nature of radiation a t lower energies has bee n 
a matter of debate JSecker et allbooi bood: IZavlin fc Pavlovl booi IKargaltsev et allbooeh 



Observations of B1929 w ith the Einst ein (Helfandlll983l), ROSAT (lYancopoulos et al. 



1994 



Becker fc Triimperlll997l ). and ASCA (jWang fc Halpernlll997l : iKawai et al.lll998l ) observato- 
ries were not able to constrain the nature of the X-ray emission from the pulsar because its 
spectrum could be equally well described as thermal or non-thermal. T he same result has 



been obtained by combining the available ROSAT and ASCA data by ISlowikowska et al. 



(J2005|), who found that the spectrum of B1929 could be fitted either by a power-law or a 
two-temperature blackbody model. 



Analyzing the images obtained in the 40 ks ROSAT PSPC observation, IWang et al. 



( I1993I ) noticed an elongated structure extendin g approximately 1 0^ behind B1929, almost 
aligned with the direction of its proper motion. I Wang et al.l (119931 ) have suggested that the 
detected tail is synchrotron radiation from the wind of B1929 c onfined by the ram pressu re. 
The detection of this diffuse emission has als o been reported by lYancopoulos et al.l (119941 ) in 
the analysis of the same ROSAT PSPC data. iKawai et al.l (|l998l ) also claim detection of some 
diffuse emission in the vicinity of B1929 in the ASCA images, but the poor angular resolution 
of ASCA prevented a detailed analysis and did not allow to separate the contribution from 
background sources in the alleged PWN emission. 



In a multi-wavelength study of B1929 (JBecker et al.ll2006l ). the tail-like structure extend- 
ing up to more than 1 pc behind the pulsar has been observed by XMM-Newton, confirming 
the ROSAT finding. The spectral analysis suggests that this extended emissio n is nonther- 



mal, l ikely produced by synchrotron radiation of the shocked pulsar wind. iBecker et al. 



(120061 ) also report on an elongated faint diffuse radio emission found in the Effelsberg ra- 
dio continuum survey data, whose brightness distribution roughly coincides with the X-ray 
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tail. These authors also conclude that the radiation of the B1929 pulsar is predominantly 
nonthermal (i.e., magnetospheric), with a rather soft power-law spectrum (photon index 
r ^ 2.7). 

To understand the nature of the extended emission behind B1929 and prove that it is 
indeed a tail of a bow-shock PWN, the nebula in the immediate vicinity of the pulsar should 
be observed with an angular resolution much better than those of ROSAT PSPC (25") or 
XMM-Newton (15" half-energy width). Therefore, we conducted two Chandra observations 
of the B1929 pulsar and its surroundings. The first results, including the de tection of a faint 



nebula in the immediate vicinity of B1929, have already been reported by iMisanovic et al. 



(120061 ). In this paper we present a more detailed analysis of these observations. 



2. Observations and results 

We observed the field around B1929 with the Advanced CCD Imaging Spectrometer 
(ACIS) aboard Chandra (Tabled]). Both observations were carried out in very faint mode 
on the ACIS-S3 chip, with the target imaged about 8" from the optical axis. The data 
were analyzed using the Chandra Interactive Analysis of Observations (CIAO) software 
(ver. 3.3.0.1; CALDB ver. 3.2.0 for observation 6657, and CALDB ver. 3.2.2 for observation 
7230). 

Table [1] also includes three archived XMM-Newton observations of B1929, which we 
analyzed in addition to the Chandra data, using SAS (ver. 7.1.0). 
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Fig. 1. — Images of the field around B1929 in the 0.3-8 keV band, combining the two ACIS- 
S3 observations. Top: Unsmoothed image (pixel size is 0.492"). The circular and elliptical 
regions mark the smallest structures detected around the pulsar at > So" levels in at least 
one of the observations (see Table [2]). Middle: The same image smoothed using a Gaussian 
of FWHM 1" . The blu e arrow shows the direction of the measured pulsar's proper motion 
( jChatterjee et al.ll2004l ). while the black arrow indicates the proper motion corrected for the 
Galactic rotation and solar peculiar velocity (see text for details). Bottom: In this image 
the pulsar contribution is subtracted (see text for details) and smoothing with a Gaussian of 
FWHM 0.5" is applied. The brightness scales in the middle and bottom panels are selected 
to emphasize different components of the extended emission in the vicinity of B1929. The 
intensity scale is linear, with a range of to 1.8 cts/pixel in the top panel, and to 0.9 
cts/pixel (middle and bottom). The position of the pulsar is indicated by a 1.5" circle in all 
panels. 
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Fig. 2. — Radial profiles of the combined pipeline- processed and subpixelized data (solid 
lines), and MARX-simulated data (dotted lines) for various blurring parameters (see text). 
The counts are measured in annular regions centered at the pulsar's position. 
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2.1. X-ray imaging 



We first produced new level -2 event files with the p ipeline pixel randomiz ation disabled 



and then applied the method by lMori et al.l (J200l|) and lTsunemi et al.l (l200ll ). which allows 



one to reach a subpixel spatial resolution. To minimize the background contribution, we 
chose the energy band of 0.3-8 keV. 

The only relatively bright source on the S3 chip is the B1929 pulsar, centered at R.A. 
= 19''32'"13.999^ decl. = 10°59'32.64", and R.A. = 19'^32°^14.000^ decl. = 10°59'32.85" 
( J2000) in the first and second observation, respectively (the centroiding 1 a uncertainty is 
0.02" for each coordinate, as determined by the CIAO procedure celldetect). The offsets 
of these positions from the radio positions of the pulsar extrapolated to the epochs of the two 
observations (0.17" and 0.14" in R.A., —0.04" and 0.16" in deck, respectively) are smaller 
than the error in the absolute Chandra astrometry, 0.45" at the 90% confidence level for 
on-axis observations on the S3 chip^. 



2.1.1. Combined ACIS image. 

Although both images show some extended emission in the immediate vicinity of the 
pulsar, its surface brightness is very low. Therefore, to increase the signal-to-noise ratio, we 
have to combine the two ACIS observations. To align the images, we reprojected the event 
files to a common position, corrected the aspect solution of the second observation for the 
pulsar's positional difference of 0.02" in R.A. and 0.21" in deck, and added the two images. 

A nebula surrounding the pulsar is clearly seen in the unsmoothed and smoothed images 
shown in the top and middle panels of Figure [H To separate this extended emission from 
the pulsar, we simulated the ACIS-S3 point source observations using the Chandra data 
simulator MARXcl. We produced simulated images of a point source at the same position on 
the detector, and with the same X-ray flux and spectral shape as those of the observed pulsar, 
running MARX for an effective exposure 100 times longer than the actual observations to 
reduce statistical errors. The simulated images were then scaled to the actual exposure times, 
combined, and compared with the data, both pipeline-processed and those with subpixel 
resolution. 



'^See §5.4 and Fig. 5.4 in the Chandra Proposers' Observatory Guide, ver. 10, at 



http://asc.harvard.edu/proposer/POGI 



See http://space.mit.edu/CXC/MARX/ 



The width of the simulated point spread function (PSF) depends on the value of the 
MARX parameter Dither Blur, which is a combination of the aspect reconstruction error, 
ACIS pixelization, and pipeline pixel randomization. According to the MARX manuao, 
appropriate values of this parameter are about 0.35" if the pixel randomization is applied 
and 0.20"-0.25" if it is switched off, but the best values of this parameter may vary from 
observation to observation. 

We found that for the default Dither Blur = 0.35", the simulated PSF is significantly 
broader than the core of the observed B1929 image, not only for the sharper image with 
subpixel resolution but also for the pipeline-processed one (see Fig. |2]). We repeated the 
simulation for a number of smaller Dither Blur values and found that the PSF core of the 
simulated data with Dither Blur reduced to 0"15 matched well to the data with subpixel 
resolution, while a simulation with Dither Blur = 0.25" is close to the pipeline-processed 
data. In both cases, the shape of the simulated PSF matched well to the observed one only 
up to ~ 0.6" from the center, indicating the contribution of extended emission at larger 
distances from the pulsar. Using the excess of the data counts with respect to the simulation 
in the 0.5"-1.5" annulus, we estimate the signal-to-noise ratio of the PWN to be > 20. For 
comparison, we also simulated the ACIS-S observation of the central compact source of the 
Cas A SNR (ObsID 6690) and found no indication of an extended emission around this point 
source. 

The PWN image in the vicinity of B1929, obtained by subtracting the simulated data 
with Dither Blur = 0.15" from the subpixel resolution image and slight smoothing with a 
Gaussian of 0.5" FWHM, is shown in the bottom panel of Figure [H The image reveals a 
compact, ~ 9" x 5", emission elongated in the direction perpendicular to the pulsar's proper 
motion, two faint wings extending in the direction opposite to the pulsar's motion and seen 
up to ~ 11" from the pulsar, and a hint of a short, ~ 3", linear structure immediately behind 
the pulsar. 



•^See also http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide/nociel6.html for a more detailed 
discussion. 
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Fig. 3. — Le/t; ACIS-S3 image of B1929 and its environment produced by combining 
observations 6657 and 7230, in the 0.3-8 keV band. The unbinned image is smoothed using 
a Gaussian of FWHM 2". The elhptical regions mark the smallest structures detected around 
B1929 (see Table [2] and Fig. [1]). The pulsar's position is marked by a 1.5" circle, while the 
40" X 20" box shows the region used for the background evaluation. The arrows indicate 
the direction of the pulsar's proper motion (as in Fig. [T]). The green circles enclose the 
optical/infrared sources in the field of view. The X marks the X-ray source found with 
CELLDETECT in the combined image, while Y and Z mark other possible X-ray point sources 
that might be just below the decection threshold. Middle: Unsmoothed combined image 
showing the same region. The image is binned by a factor of two. The ellipse shows the 
region of enhanced emission, in which the diffuse emission is detected at a ~10cr level in 
the combined image. The red box and the line mark the size of the rectangular region and 
the direction in which the region is moved to produce the linear profile shown in Fig. |H 
Right: The same unbinned image smoothed using an adaptive Gaussian kernel with the size 
self-adjusting to show the structures with the signal-to-noise ratio in the range 2.2 to 4. 
The overlaid contours correspond to the combined XMM-Newton observations. The contour 
levels are the same as in Fig. O The intensity scale is linear in all images, with a range of 
to 0.2 cts/pixel in the left and right panels and to 1.8 cts/pixel (middle). 
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After confirming that there is a PWN in the immediate vicinity of B1929, we proceeded 
by searching for the nebular emission farther from the pulsar. Figure [3] shows a 2.5' x 2.5' 
image combining the two Chandra observations. The image on the left is smoothed using a 
fixed-size Gaussian of 2" FWHM, while the right panel shows a heavily smoothed image using 
an adaptive kernel with an adjustable size to show the structures with the signal-to-noise 
ratio in the range 2.2 to 4. The middle panel shows the combined unsmoothed data binned 
by a factor of two. The images show several regions of extended emission: a bow-shock- 
shaped extended emission with the Southern Wing (region 1) and Northern Wing (region 
2), some faint emission in front of the pulsar (region 3 or the Front), faint diffuse emission 
elongated in the direction opposite to the direction of the pulsar's projected velocity, with 
an enhancement at the distance of ~30" (region 4 or the Inner Blob). A more extended but 
fainter region 5 (Outer Blob) is seen farther behind the pulsar, at a distance of approximately 
1.5'-2'. Regions 1 through 3 are also shown in the top panel of Figure [H while the numbers 
of counts in the regions are given in Table [2l 

Although these structures (regions 1 to 5) were detected at a level of at least 3a in one 
or both individual observations, the small number of detected counts prevented the detailed 
spectral analysis. Therefore, we selected a larger region of the extended emission to study 
the spectral properties of the PWN ( §2.3.21) . The region is marked by the ellipse ( semi-axes 
12" and 25") in the middle panel of Figure [3l detected at a level of at least 5a in each 
observation. 

The blue arrows in Figures [H and [3] indicate the direction of the pulsar 's proper motion 



(fin = 94.09 ± 0.11 mas yr ^, fis = 42.99 ± 0.16 mas yr ^) as measured by IChatterjee et al. 



(I2OO4J ) with respect to extragalactic reference radio sources. To determine the proper motion 
with respect to the local ISM, we corrected the measured proper motion for the effects of 
Galaxy rotation and the Sun's pec uliar motion with respect to the local standard of rest 
(LSR). Following the procedure by 
for the Sun's peculi a r velo city from 



Johnson fc Soderblom 



Il987h . and adopting the value 



Boesgaard fc Tripiccd ( 19861 ) and Oort constants from 



Binney fc Tremaind (119871 ). we obtained the corrected proper motion Ha = 91.6 ± 0.2 mas 
yr~^, ^s = 49.1±0.2 mas yr~^, shown by the black arrows in Figures [1] and [3l The total proper 
motion, n = 104.0±0.3 mas yr~^, corresponds to the transverse velocity v± = 178±8 km s~^. 
Its direction (position angle 61.7° ± 0.1°, counted East of North) seems to be better aligned 
with the extended tail of diffuse emission, which further supports the PWN interpretation 
of the detected X-ray emission around this pulsar. 

Since the observed diffuse emission might be partly due to unresolved point X-ray sources 
in the field of view (e.g., stars or background AGNs), we examined available optical, near- 
infrared, and radio data of the region around B1929. We found several sources (marked with 
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the g reen circles in Fig.[3l left) in the USN0-B2 flMonet et al.ll2003r ) and 2MASS flCutri et al. 



20031 ) catalogs, which coincide with the regions resembling faint point-like objects in our 
smoothed images. Although these possible sources are very faint in the Chandra images 
(only one of these sources, marked with an X, was detected by the CIAO tool celldetect 
above the source threshold of 3, while Y and Z mark other two possible X-ray point sources 
that might be just below the detection threshold), they might contribute to the observed 
extended emission. It is clear, however, that most of the detected tail-like structure is indeed 
diffuse emission associated with the pulsar. 

Figure H] shows the brightness distribution along the tail, which further demonstrates 
the excess emission over the background in front of the pulsar, and also in the direction 
approximately opposite to the projected pulsar's velocity. The photons were collected from 
1" X 35" rectangular regions (see Fig. E]) along the tail, and also from the background regions 
of the same area on both sides of the tail in the unsmoothed combined Chandra image. We 
show only the brightest part of the tail, up to ~ 40" from the pulsar, because the extended 
emission farther away is too faint and cannot be distinguished from the background using 
such small collecting regions, but it is much better seen in the smoothed image. 
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Fig. 4. — X-ray linear brightness profile along the direction of the pulsar's tail. The black 
data points were obtained by counting the photons in the 1" x 35" rectangular regions (see 
Fig. [3]). The pulsar moves towards the negative values on the X-axis. The red data points 
show the X-ray profile of the background emission, evaluated in the rectangular regions of 
the same size on both sides of the tail emission. 
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To look for extended emission at even larger dist ances from the pulsar and compare it 
with the 10' — 15' tail detected by lBecker et al.l (120061 ) . we produced heavily binned images 
shown in Figure O Although the smaller ACIS-S field of view and the relatively short 
Chandra exposure do not allow us to detect the full extent of the faint diffuse emission 
observed by XMM-Newton, at least some parts of the B1929 long tail are seen up to ~ 4' 
from the pulsar in the combined ACIS-S3 image. The long tail extends from the compact 
PWN resolved in the Chandra images in the same direction as the emission detected in 
XMM-Newton data. 
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Fig. 5. — 10' X 8' ACIS-S3 field around B1929. The image on the left is produced from the 
combined event file, binned by a factor of 4 and then smoothed by a Gaussian of FWHM 
6". The contours, at levels 0.08, 0.1, 0.12, 0.13, 0.14, 0.15, 0.2, 0.3, 0.5, 0.9, 2 and 5 
counts pixel"^, correspond to the three combined XMM-Newton observations (MOSl and 
M0S2 data, pixel size 1.1"). The right panel shows the same image, which is binned by a 
factor of 8 but not smoothed. To minimize the background, both images include only the 
0.5-7 keV energy range. The intensity scale is linear, with ranges of to 0.8 counts pixel"^ 
(left) and to 6 counts pixel"^ (right). 
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2.1.2. Individual observations and variability 
Since some PWNe sho w significant tempo r al variations on tim e-scales of months (e.g 



the Crab and Vela PWNe, iHester et al.l l2002l : iPavlov et al.ll2003l ). we have examined the 



individual ACIS-S observations, which are separated by ~6 months, in search for variability 
of the PWN of B1929. Figure [H] shows the images of the individual Chandra observations of 
B1929 and its PWN. Similarly to the combined data, the images on the left were smoothed 
using a Gaussian of FWHM 2", while on the right we show the same images adaptively 
smoothed to reveal structures with the signal-to-noise ratio in the range 2.2 to 4. The bright 
features seen in the combined image (Fig. [3]) are also visible in these individual images, but 
their apparent brightness seem to show a substantial difference between the observations. 
For example, the Inner Blob (region 4), clearly detected in the observation 6657 (December 
2005), is not visible in the image of the observation 7230 (May 2006), while the Outer Blob 
(region 5) appears fainter and slightly smaller in the second observation. In addition, the 
Southern Wing and the emission in front of the pulsar (region 3) are slightly brighter in the 
observation 6657, while the Northern Wing is brighter in the second observation. 

To quantify this apparent variability, we measured the total background-subtracted 
counts detected in each of the extended regions in the unsmoothed images of both data sets, 
and calculated the corresponding surface brightness (listed in Table [2]). The background 
emission was determined using the rectangular region with an area of 784 arcsec^ shown in 
Figure [3] (left), from which we measured (3.3 ± 0.4) x 10~^ ctss~^ arcsec"^ and (2.2 ± 0.4) x 
10~^ctss~^ arcsec"^ in the observations 6657 and 7230, respectively. 

As shown in Table [2|, the surface brightness of the five extended regions is found to 
change by a factor of two (e.g. region 1) or even four (region 3). However, due to the large 
statistical errors, the variability could only be detected at a level of up to ~ 2.5cr (region 5). 
Deeper observations would be needed to further examine possible variability of the nebula 
associated with this old pulsar. We also note that the variation of the pulsar count rates 
between the two Chandra observations was found to be within the statistical uncertainty 
(see § [MID . 
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Fig. 6. — Le/t: 0.3-8 keV images of the first (top) and second {bottom) ACIS-S3 observations, 
smoothed using a Gaussian of FWHM 2". Right: The same images adaptively smoothed to 
show the structures with the signal-to-noise ratio in the range 2.2 to 4. The intensity scale 
is linear with a range of to 0.2 cts/pixel in all panels. 
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2.2. Connection with the XMM-Newton data 

To examine the connection between the Chandra and XMM-Newton data on the PWN, 
we produced combined MOSl and M0S2cl broad-band images for eac h of the three avail able 



XMM-Newton observations of the B1929 field listed in Tabled] (see iBecker et al.l l2006l . for 
more detail about the XMM-Newton data and for the combined image of all three observa- 
tions). 

The smoothed images are shown in Figure [71 overlayed with the structures detected 
in the ACIS-S observations. Since the MOS PSF is much broader than that of ACIS, the 
compact PWN (regions 1-3 in Figs. [1] and [3]) and the beginning of the PWN tail (region 
4) cannot be resolved from the pulsar. However, at least in the first two XMM-Newton 
observations we clearly see the local brightening in the PWN tail ~ 2' from the pulsar, 
approximately coinciding with the Outer Blob (region 5) detected in the Chandra data. 
From the three unsmoothed MOSl images, we calculated an average background-subtracted 
count rate of 1.5 ±0.3 counts ks~^ in the Outer Blob. With account for different ACIS-S and 
MOSl responses, this count rate is consistent with that measured in the combined Chandra 
data in the same region (see Table [2]). 



*The PN observations were taken in the small-window mode and have high background. 
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Fig. 7. — Combined M0SH-M0S2 images of B1929 and its surroundings in the three XMM- 
Newton observations. Broad-band (0.3-10 keV) images were smoothed using a Gaussian of 
FWMH 6". The five regions of enhanced emission and possible pointhke sources X, Y and 
Z detected in the ACIS-S observations are shown (see left panel of Fig. [3] and Table [2]). 
The bright source (Var) ~ 2.5' south from the pulsar exhibits strong variability, also clearly 
visible in the unsmoothed MOS images. The intensity scale is linear, with a range of to 6 
counts per pixel in all panels. 
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Since our Chandra data indicate possible variability of the PWN, we examined the 
XMM-Newton data for brightness variations in the extended emission. From the number 
of the background-subtracted counts in the Outer Blob, we determined the averaged MOSl 
count rates of 1.7 ± 0.5, 1.6 ± 0.5, and 1.3 ± 0.5 counts ks~^ in observations 0113051301, 
0113051401 and 0113051501, respectively. Thus, although the Outer Blob changes its surface 
brightness by a factor of 2 at a 2.5cr level in our ACIS-S data, the surface brightness of the 
same region in the MOSl data, separated by approximately the same period of ~5 months 
as the two Chandra pointings, does not show a significant brightness change. 

Figure [7] suggests that the region southwest of the Outer Blob might be variable. How- 
ever, in the ACIS image we see no diffuse emission in this area, but there are two faint 
pointlike sources, X and Y (see Figs. [3] and [7]). Hence, the apparent variability of this region 
could be due to the variability of these two sources, which could not be resolved in the XMM- 
Newton observations. Thus, although the Chandra observations hint that the tail of B1929 
might be variable on time-scales of several months, longer high-reslolution observations are 
required to firmly establish this variability. 

Interestingly, the XMM-Newton images show a highly variable source (Var) ~ 2.5' south 
of the pulsar. It was almost as bright as the pulsar in the 2004 April 27 observation, while 
it is barely visible in the 2003 November 10 image. The source was at the edge of the 
S3 chip in the first Chandra observation and barely detectable. However, it completely 
disappeared in the second observation, although its position was within the field of view. 
As it is apparently a point source, strongly detached from the other PWN structures, we 
conclude that it is unrelated to the B1929 PWN. 
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Fig. 8. — Top: Spectra of B1929 extracted from the Chandra observations 6657 (black) and 
7230 (red) and the best-fit absorbed PL model. The model parameters are given in Table [31 
Middle: PN, MOSl and M0S2 spectra of B1929 extracted from the three XMM-Newton 
observations (PN: black, red, green; MOSl: blue, cyan, magenta; M0S2: yellow, orange, 
lime-green; in the observations 0113051301, 0113051401 and 0113051501, respectively), and 
fitted simultaneously with the Chandra data (Obs. 6657: dark-green; Obs. 7230: dark- 
blue) using the absorbed PL+BB model. The model parameters are given in Table [31 
Bottom: Photon spectrum for the best-fit absorbed BB+PL model and its components. The 
unabsorbed luminosity of the BB component, ~ 1.2 x 10^° ergs s~^, is estimated to be ~45% 
of the total luminosity, ~ 2.6 x lO^'' ergs s~^, in the 0.3-8 keV band. 
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Fig. 9. — Top: 68%, 90%, and 99% confidence contours for tlie Chandra ACTS and XMM- 
Newton EPIC spectra of B1929 (dashed and solid lines, respectively), computed for two 
interesting parameters for the PL model. Bottom: 68%, 90%, and 99% confidence contours 
for the Chandra and XMM-Newton observations fitted simultaneously with the PL model. 
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Fig. 10. — 68%, 90%, and 99% confidence contours for various parameter pairs for tlie 
PL+BB fit to the spectrum of B1929. The confidence levels correspond to two interesting 
parameters. The Chandra ACIS-S and XMM-Newton spectra are fitted simultaneously with 
all parameters allowed to vary. The top panel shows the lines of constant bolometric lumi- 
nosity of an equivalent sphere (see § 3.2) in units of 10^° ergs s^^ for the BB component, while 
the middle panel shows the lines of constant unabsorbed fiux in units of 10~^^ergscm^^ s~^ 
for the PL component. 
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2.3. Spectral analysis 

2.3.1. The pulsar 

The spectral Chandra data reduction was done by applying the standard procedures in 
CIAO. From a circular region with a radius of 1.6" (or 3.3 pixels) centered on the pulsar, we 
extracted 645 and 690 counts from the first and second observation, respectively. Four circu- 
lar background regions, each with a radius of 5" , were selected in the vicinity of the pulsar 
and then combined and used for the background subtraction. The background contribution 
in the source aperture is negligibly small (~ 1 count in each of the observations), but the 
aperture may include a small number of photons emerging from the compact nebula. Using 
the radial profiles in Figure [21 we estimated the PWN contribution to be around 2%-3%; 
hence it does not influence the results of the spectral fitting significantly. We grouped the 
extracted counts in a minimum of 15 and 17 counts per energy bin for the spectral analysis. 
The response matrices (rmf) were calculated using the CIAO task mkacisrmf, and the 
spectral analysis was performed using XSPEC (ver. 11.3.2). 

Since we did not find statistically significant changes between the two Chandra obser- 
vations, we fit both data sets simultaneously. The absorbed power-law (PL) model yielded 
a very good fit (Tabled Fig. [8], top), albeit with surprisingly large best-fit hydrogen col- 
umn density, A^h ~ 2.5 x 10^^ cm~^, and photon index, F si 2.9. On the contrary, we were 
not able to produce a statistically acceptable fit with the absorbed blackbody (BB) model 
as the discrepancy between the best-fit model and the data was large below 0.7 keV and 
particularly above 2 keV, resulting in x^ = 1.86 (vs. x^ = 0.91 for the PL fit) for 71 d.o.f. 

The pulsar models predict that a thermal component from heated polar caps, and 
perhaps from the whole NS surface, may be detectable, making the main contribution at 
lower photon energies. Therefore, we attempted to fit the Chandra spectra with the PL-fBB 
model. The addition of the BB component, with the best-fit temperature kT Ri 0.3 keV and 
emitting area A ~ 4 x 10^ m^, resulted in a substantially smaller hydrogen column density. 
Ah < 1.4 X 10^^ cm~^, and a softer PL spectrum, F ^ 1.8 (see Table [3]), but the probability 
that the BB component is required by the data was only 69.5%, according to the F-test. 
We also tested the BB+BB model, but it yielded either statistically unacceptable fits or 
unphysical fitting parameters (e.g., kT ~ 100 keV, or an extremely low absorption column). 

To compare the Chandra ACIS spectra with those obtained with the XMM-Newton 
EPIC detector, we extracted and fitted the pulsar spectra from the three XMM-Newton 
observations. The source counts were extracted from a circular region with a radius of 
13" centered on the pulsar, while a 38"-radius circular region northwest of the pulsar was 
used for the background subtraction. We collected a total of 450 (with 96.8% of the source 
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contribution) , 658 (94.8%) and 805 (75.3%) counts in the PN (small- window mode) observa- 
tions 0113051301, 0113051401, and 0113051501, respectively. In addition, we extracted 172 
(98.0%), 282 (97.3%), and 110 counts from the MOSl, and 153 (97.3%), 276 (98.0%), and 
121 (96.8%) from the M0S2 observations, with MOS cameras operated in the large window 
mode. The EPIC response matrices and ancillary files were calculated using the SAS tasks 
RMFGEN and ARFGEN. From the number of the PWN counts detected with ACIS in an 
aperture of 13", we estimate that the contribution from the nebula in the EPIC spectra is 
about 5%. 

Fitting all the 9 XMM-Newton data sets (PN, MOSl and M0S2 in the three observa- 
tions) simultaneously, we found that the PL model produced a statistically acceptable fit, 
with the model parameters consistent with those obtained from the Chandra data (Fig. ^^ 
Adding the BB component to the spectral model resulted in about the same fitting parame- 
ters as for the Chandra data alone (Table [3]), but improved the fit more significantly (94.8%, 
according to the F-test). 

Finally, we fit the Chandra and XMM-Newton data simultaneously, first with the PL 
model. Because of the larger PWN contribution in the XMM-Newton data, we first linked 
all the model parameters except for the PL normalization. However, since we found that the 
difference between the PL normalizations in the Chandra and XMM-Newton spectra was 
within the statistical errors, we finally linked all the model parameters. As expected, we 
obtained the PL parameters similar to those obtained from the separate fits to the Chandra 
and XMM-Newton data, but with smaller uncertainties (Fig. [HI bottom). 

Fitting the combined data with the PL+BB model, we found that the addition of the 
BB component provided a very significant improvement of the fit (99.98%, according to the 
F-test). The best PL+BB fit is shown in Figure El the fitting parameters are given in Table [31 
and the confindence contours for various model parameters are presented in Figure dUl We 
note that the hydrogen column density, A^h < 4 x 10^°cm~^, is significantly lower in the 
PL+BB fit and consistent with the value estimated from the pulsar's dispersion measure, 
while the slope of the PL components, P ~ 1.7, is similar to those found in many other 
pulsars (see § 13.21 for further discussion). 



^We should mention that our PL fitting parameters differ significantly from those obtained by Becker 
et al. (2006) from fitting the combined XMM-Newton EPIC and ROSAT PSPC data (they do not provide 
model parameters from fitting the XMM-Newton data separately). For instance, these authors obtained a 
lower hydrogen column density, A'h = (1-6 ± 0.2) x 10^^ cm~^, and a smaller photon index, F = 2.72J;Qj!,g. 
We believe that the discrepancy is caused by systematic errors due to poor cross-calibration of the PSPC 
and EPIC detectors. 
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Fig. 11. — Toj): 68%, 90%, and 99% confidence contours in the A^h — T plane of tlie PL 
model for tfie PWN spectral data (the elliptical region in Fig. [3]). Bottom: 68%, 90%, and 
99% confidence contours in the F-Normalization plane of the PL model with the absorption 
column fixed to the best-fit value obtained from the spectrum of the pulsar. 
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2.3.2. The nebula 

As the detected PWN is very faint, we have to choose a sufficiently large extraction 
region to measure its spectrum. The large-scale PWN appears to consist of a smaller, 
slightly brighter part up to 30" behind the pulsar (marked by the ellipse in the middle panel 
of Figure [3]), and a fainter emission farther down, up to a distance of ~2'. From the brighter 
part of the nebula region (the ellipse) we collected 138 and 136 counts (64% and 63% of 
which are from the nebula) in the ffist and second observation, respectively. 

To constrain the model parameters, we fit the spectra from the two observations simul- 
taneously, with the absorbed PL model. The fit yielded a PL slope r=1.0-1.7 at a 90% 
confidence level (Fig. [11], top) and an absorption column N^ = 0.0 — 0.9 x 10^^ cm~^ (90% 
confidence). We also fitted the PWN spectra with the absorption column fixed to the best-fit 
value of A^H = 1-7 x 10^^ cm~^ obtained from the spectrum of the pulsar, but the photon 
index range remains the same (Fig. [TT], bottom). We also attempted to fit the spectra from 
the PWN region farther down the tail, to see if there is any indication of spectral softening, 
but the number of photons from this faint part of the nebula is too small for a miningfuU 
spectral analysis. 

The PWN flux, -FpwN = 2.1-3.6 x 10^^^, determined from the elliptical region, is poorly 
constrained because of the low statistics. We also measured a flux of 3.4 — 5.8 x 10"^*^ 
ergs cm~^ s~^ from the whole visible PWN. The unabsorbed luminosity of the whole nebula 
estimated from this observed flux is LpwN = 5.3-8.6 x 10^^ ergs s~^, in the 0.3— 8keV band. 

The absorption column determined from the spectral fit of the nebula is significantly 
lower than that determined from the PL fit of the pulsar spectrum, but it is in a good agree- 
ment with the value derived from the PL-I-BB model. This further supports our detection 
of the thermal component in the spectrum of B1929. 



3. Discussion and conclusions 
3.1. The nebula around PSRB1929+10 

3.1.1. P WN morphology 

In the analysis of the two ACIS-S observations of B1929, we detected a faint PWN 
surrounding this old pulsar. We found that the pulsar is immersed in a very compact, 
K. 9" X 5" nebula elongated in the direction perpendicular to the pulsar's proper motion 
(Fig. 1, bottom). Together with the two patchy "wings" (best seen in slightly smoothed 
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images), the overall compact structure looks like a bow shock PWN, with an apex at ^ 3" 
[^ 1.6 X 10^^ cm at c? = 361 pc) ahead of the pulsar. In addition, there is marginal evidence of 
a short jet-like feature emerging from the pulsar between the wings, in the direction opposite 
to the pulsar's proper motion (seen up to 3"-4" from the pulsar in the unsmoothed combined 
image and in the pulsar-subtracted image in Fig. [1]). There is also some very faint emission 
in front of the apparent bow-shock apex, seen up to 6"-7" from the pulsar (see Fig. Hj). 

In addition to the compact PWN, a faint, inhomogeneous tail-like structure was de- 
tected, extending in the direction approximately opposite to the pulsar's proper motion up 
to at least ~ 1.5'-2'. This structure is well aligne d with the rnuch longer tail (up to 10^ -150 
detected in the ROSAT and XMM-Newton data JWang et al.lll993l : ISecker et allbooeh . 



The surface brightness of the detected PWN is not uniform: e.g., the two wings are 
brighter than the other PWN structures (see Table [2]). Furthermore, the Chandra and XMM- 
Newton data suggest a possible variability of the detected PWN, although the number of 
detected photons was too small to establish it unambiguously. 

The approximate axial symmetry with respect to the direction of the pulsar's proper 
motion strongly suggests that the PWN is shaped by the oncoming flow of the ambient 
matter in the pulsar's reference frame. The corrected pulsar's proper motion and distance 
correspond to the pulsar's speed fpsR = 178 (sin i)~^ km s~^, where i is the angle between the 
velocity vector and the line of sight. It exceeds the ISM sound speed, Cg = {5kT/3fimHY^'^ = 
STfi^^^'^T^ km s~^ (where fx is the molecular weight, T5 = T/10^ K), for T < 2.3 x 
10^/x (sini)~^ K, i.e., for virtually any plausible ISM temperature, which supports the bow- 
shock interpretation of the detected nebular emission. 

Our current understanding of the pulsar wind dynamics in a bow-shock PWN mostly re- 
l ies up o n numerical models. A r ecent work in this dir e ction has been done by Ivan der Swaluw 
(l2003l ). iBucciantini et al.l (120051 ) . and IVigelius et al.l (120071 ). These models show that at su- 
personic pulsar speeds the termination shock (TS) of an isotropic pulsar wind acquires a 
bullet-like shape, with a distance 



Rh ~ {E/Aticp, 
between the pulsar and the bullet head, where 



a/2 



Pram — Pamb ^PSR 



(2) 



is the ram pressure, and pamb is the ambient density. The shocked pulsar wind is confined 
between the TS and the contact discontinuity (CD) surface, while the forward bow shock 
(FBS) separates the shocked ambient medium (between the CD and the FBS) from the 
unshocked one. For large Mach numbers, Ai = Vpsn/cg = (3pi.am/5pamb)^^^ ^ 1, and small 
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values of the magnetization parameteio of the pre-shock pulsar wind, a < 0.1, the bullet's 
cylindrical radiu s is Tts ~ Rh and the dis tance of its back surface from the pulsar is R^ ~ QRh 



(see Figs. 1-3 in iBucciantini et al.ll2005l ). The shape of the CD surface ahead of the pulsar 
is similar to that of the TS but with the apex at Rqt) ~ 1-3-Rh; while it acquires a cylindrical 
shape with a radius rcD ~ ^Rh behind the TS bullet. The shocked wind flows away from 
the pulsar within this cylinder, forming a PWN tail. The flow velocity in the central part 
of the tail (the inner channel, r < rxs) is ~ 0.1c-0.2c, while it is as high as ftaii = 0.8c-0.9c 
in the bulk of the tail's volume, at tts ^ '^ < '"cd (where r is the cylindrical radius). The 
magnetic field, purely toroidal in these models, also depends on r, being enhanced toward 
the tail axis and the CD surface. Its typical value in the tail is 

5taii ~ 2.8 (c/t;taii) (ap,am)'/'G^ (3) 

where the pram is in units of dyn cm~^ (see eq. (13) of IBucciantini et al.l l2005l ). In the 



numerical simulations, the flow velocity and the magnetic field do not show substantial 
changes along the tail, but these simulations are l imited to relatively short distances from 



the pulsar (e.g., < 13i?/j in IBucciantini et al.ll2005l ). 



To interpret the B1929 PWN in the framework of the current MHD models, we at- 
tempted to match the observed morphological structures with those predicted by the sim- 
ulations. As the brightest X-ray emission is expected from the shocked wind between the 
TS and CD in the front part of the PWN, we assume that the very compact structure in 
the immediate vicinity of the pulsar (Fig. (H bottom) is the PWN head and estimate the 
distance from the pulsar to the CD apex to be ~ 3". This corresponds to Rct> ~ 1-6 x 10^^ 
cm, Rh ~ rTs ~ 1-2 x 10^*^ cm (or 2.3"), Rb ~ 7.5 x 10^^ cm (or 14"), and rco ~ 4.8 x 10^*^ 
cm (or 9"). In this picture, the patchy wings in the PWN image are produced by the 
shocked wind between the TS and CD, while the "blobs" farther behind the pulsar are in 
the pulsar tail. However, the faint emission ahead of the pulsar (region 3, or the Front, in 
Fig. [1], top, and Table [2]) is not explained in this model. We can only speculate that this 
emission might be produced by a polar outflow along the pulsar's spin axis (not included in 
the models), possibly scattered by the head wind, or it might suggest that, because of the 
Rayleigh- Taylor instability, the CD surface acquires an irregular shape ahead of the pulsar, 
with "flngers" of the shocked wind penetrating into the shocked ISM. As the above estimates 
of sizes of the PWN elements are based on simplifled models and shallow images, one should 
consider them as order-of-magnitude estimates rather than accurately established values. 
Therefore, in the following estimates, we will retain the explicit dependence on i?/i, scaling 



it as Rh = Rh M^ x 10^^ cm. 



^The magnetization parameter in the pre-shock wind is defined as the ratio of the Poynting flux to the 
kinetic energy flux. 
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Using this estimate for Rh and equation (1), we can calculate the ram pressure: 

Pram = E fn{^T:cRl)~^ = 1.0 X 10'^%^M6 ^JH cm'^, (4) 

where we introduced the factor /j^ to account for possible anisotropy of the pulsar outflow 
(e.g., /q < 1 for a mostly equatorial outflow if the equatorial plane is perpendicular to the 
pulsar's velocity). Since, on the other hand. 

Pram = PambV^PSR = 5-3 X 10~^°namb (siui)"^ dyu cm~^ , (5) 

where riamb = Pamb/'"^H, we obtain an estimate for the ambient number density: 

namb = Q.l'^fnRh% sin^ i cm"^ . (6) 

Using equation (4), we can estimate the Mach number: 

M = 7.7/^/^i?,->Xi2 , (7) 

where we scaled Pamb to 10^^^ dyn cm^^, a typical value for the ISM. From the corresponding 
sound speed, Cg = '23 f^ -R/i.iePamb -i2(si'^^)~^ ^^ ^~^ ^^ estimate the temperature: 

T = 3.9 X lOV/n'^' i6Pamb,~i2(sin2)-2 K . (8) 



The arnbient pressure in the unperturbed ISM is in the range Pamb,-i2 ~ 0.2-2 (JHeilesll200ll : 



Ferrierell200ll ). The pressure in the pulsar's neighborhood might be somewhat higher if the 



pulsar's UV and soft X-ray emission heats and ionizes the surrounding medi um, which is 



suppo rted by the lack of an Hp. bow-shock PWN in the observations reported by lBecker et al. 



(120061) (see iBucciantini &: Bandieral l200ll . for a detailed discussion on conditions of observ- 
ability of Ha bow-shock PWNe). As the factor f^f^^R^ ]^g(sinz)~^ is also very likely greater 
than unity, we expect T ~ 10^ K to be a realistic estimate for the ambient temperature. 

We should note that although the comparison of our data with the numerical mod- 
els leads to reasonable estimates for the parameters of the ambient medium, these two- 
dimensional models are based on a number of assumptions that are not necessarily realistic. 
In particular, the models assume an isotropic pulsar outflow, whereas we know that it is 
predominantly equatorial, at least in the case of young pulsars. Obviously, a substantial 
anisotropy would distort the shape of the bow shock in the pulsar vicinity (e.g., decreasing 
Rh in the case of an equatorial outflow in the plane perpendicular to the pulsar velocity). 
For B1929, a hint of an equatorial outflow might be suggested by the elongation of the most 
compact PWN structure perpendicular to the proper motion (Fig. (H bottom). Such a struc- 
ture, together with the "wings" , might be interpreted as a torus of the shocked equatorial 
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wind, compressed and possibly bent by the oncoming wind. Such an interpretation might 
also be supported by the marginal detection of the jet-like structure immediately behind the 
pulsar and the faint emission ahead of the pulsar (a crushed counter-jet?). 

Application of the PWN simulations to the observed tail leads to some discrepancies. 
For instance, although the tail may look like a cylinder with the predicted diameter of ~ 18" 
at the relatively small distances, < 2', from the pulsar in the ACIS images, its transverse size 
apparently increases at larger distances, possibly reaching ~ 5' at the edge of the MOS field 
of view (15' from the pulsar). The tail's surface brightness looks very patchy, which might 
suggest some interaction of the shocked wind with the surrounding medium, such as the 
Kelvin-Helmholtz instability caused by the shear at the CD surface. This instability may 
produce clumps of shocked ISM embedded i n the shocked wind, wh ich would distort the 



flow pattern and compress the wind material (JBucciantini et al.ll2005l ). The same instability 



may result in temporal variations in the X-ray emission of the tail, with a time scale of 
^cd/c ~ 0.5 months, which might explain the alleged variability of the PWN (§ 12.1.21) . 

Using equations (3) and (4), the magnetic field in the tail can be estimated as 

Staii = hQa^'\Q.hc/v,^^)fll^Rl\^ fiG . (9) 

We should note that this field becomes rather small, comparable to the interstellar magnetic 
field, at the usually assumed small values of the magnetization parameter, a < 10^^. 

Using this estimate of -Btaib we can estimate the synchrotron cooling time for relativistic 
electrons in the tail, Tgyn = 5.1 x lO^^'j^^ BZ^ s, where 7 = 10*^78 is the electron Lorentz factor 
and i?_5 = B/10^^ G. As the characteristic energy of synchrotron photons is 

E ~ 0.5B_5-fl keV, (10) 

photons with energy E are emitted by electrons with 73 ~ 2(i5/lkeV)^/^-B„5 , and the 
synchrotron cooling time corresponding to E = 8 keV photons (upper energy of our band) 
is 

Tsyn = 6.8 X 10\afn)-'/\v^o./0.5c)'/\E/8keV)-'/'Rl% s. (11) 

It follows from this equation that the projected tail length, /tail ~ ^flowTsyn sini, as observed 
at energy E, can be estimated as 

'••"~i"i° [oTc) (-^j [-7m) "'"'=""■ <''> 

— 1/2 

For instance, for a = 0.01, Rh,i6fn = 2, i'aow = 0.5c, sini = 0.5, and E < 8 keV we obtain 
^taii ^ 150 pc (i.e., ~ 24°), two orders of magnitude larger than observed by XMM-Newton 
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and ROSAT. The main reason for this discrepancy is the assump tion that the flow speed 
remains very high along the entire tail in the ideal MHD model. iBucciantini et al.l (120051 ) 
have noticed that a similar discrepancy arises when the same model is applied to the X-ray 
tail of the Mouse PWN: the observed tail length is much shorter than that estimated from 
the model at any a < 1. These authors suggest that the flow can be slowed down by the 
interaction with the ambient medium, which could also explain the divergence of the Mouse 
tail in the radio. As these factors have not been included in the current models, we can only 
empirically estimate an average flow velocity that would be consistent with the observed 
length of 1.5 pc: fflow ~ 0.1ccr_2 {fn/R\ i6)^''"'^°(sini)~^/^, assuming that equation (3) is still 
applicable (here cr_2 = o"/10~^). We note that this velocity is still much larger than the 
pulsar velocity. This means that the equation 



^tail — "WPSR.XTf 



syn 



0.3(t;psR,x/178 km s-^)%^Bzl pc. 



(13) 



often used for estimating the magnetic field (e.g., iBecker et al.ll2006l ). is inapplicable in this 



cas< 



a 



It is worthwhile to mention that the length of the tail of B1929 measured so far might 
be limited by the field of view of the available XMM-Newton observations, and further deep 
observations along the tail are needed to determine its full extent. Recent observations with 
Chandra and XMM-Newto n have discovered a number of other PWNe with very long tails 



flKargaltsev fc Pavlovll2008h . includin g the 6 pc tail (k.n IRk ~ 600) behind PSR J1509-5850, 



the longest pulsar tail known to date (JKargaltsev et al.ll2008l ). Detection of more such objects 



would help to facilitate further modeling that would account for cooling in the varying 
magnetic field at larger distances from the pulsar. 



3.1.2. X-ray luminosity and spectrum of the PWN 

From our estimate of the total unabsorbed flux of the B1929 PWN detected with Chan- 
dra, we estimated the X-ray PWN luminosity, LpwN = (5.3-8.6) x 10^^ ergs s^^, which cor- 
responds to the PWN efficiency //pwn = -^pwn/-E' = (1.4-2.2) x 10""^ in the 0.3-8 keV band. 
Similar values of the luminosity and efficiency have been measured in the XMM-Newton 



^Equation (13) implies that the shocked wind flow acquires the speed of the ambient ISM matter in the 
immediate vicinity of the pulsar, so that the observed elongated PWN is actually a "trail" of the decelerated 
wind left beh ind the moving; pulsar. Such an assumption str o ngly contradic t s all t he models of bow-shock 
PWNe (e.g., iRomanova et al.1 l2005t IBucciantini et al.l 120051 ) . iBecker et al.l (|2006l) obtained a reasonable 
estimate for the magnetic field, i? < 12 /xG, using eq. (13) because they assumed 7 = 10^. However, at such 
B and 7, the energy of a synchrotron photon is _E = 6 eV (see eq. [10]), well below the X-ray band. 
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data by iBecker et al.l (120061 ) from a larger area farther southwest from the pulsar. As both 
measurements were taken for small fractions of the PWN, the luminosity and efficiency of 
the entire PWN can be higher than these estimates. 

The estimated efficieri c y of t he B1929 PWN is within the range t^pwn ~ lO^^'^-lO^^ 
found by iKargaltsev et al.l (120071 ) in a recent study of several energetic middle-aged Vela- 
like pulsars and their PWNe. Compared to a few known PWNe with long tails, which, 
on average, show high er efficiencies, '/^pwN ~ lO^^'^-lO^^'^, than more compact PWNe (see 
Kargaltsev et al.ll2008l ). the B1929 PWN is among the least efficient ones; however, we stress 
that a deeper observation of B1929 is needed to measure the efficiency more accurately. 

The spectral slope o f the PWN spectrum is also consistent with the v alues measured 
for other PWNe (e.g., see IKargaltsev et al.l 120071 : IKargaltsev fc Pavlovll2008l ). 



3.2. Spectral properties of PSRB1929+10 

As we mention in ^ it has been a matter of debate whether the X-ray emission from 
B1929 is predominantly magnetospheric or it has a signifi c ant th ermal component emitted 
from the neutron star surface. Most recently, iBecker et al.l (|2006l ) have concluded, based on 
the analysis of the XMM-Newton data, that the spectrum of B1929 is best described by a 
PL model (i.e. the emission is predominantly magnetospheric), while the contribution of the 
thermal component, modeled as BB radiation, does not exceed ~7%. Our analysis of the 
Chandra and XMM-Newton spectra has shown, however, that adding the BB component to 
the PL model significantly improves the fit to the combined Chandra plus XMM-Newton 
spectrum, with the the best-fit BB component providing ~45% of the luminosity in the 
0.3-10 keV band. 

In addition to the better fit, there are other serious arguments in favor of the PL+BB 
model. First of all, the hydrogen column density, A^h = (2.23 ± 0.27) x 10^^ cm~^, obtained 
from the PL fit, looks unreasonably large. In particular, it is much larger than the standard 
estimate, A^h ~ 10A"e ~ 1 x 10^° cm~^, obtained from the pulsar's dispersion measure, DM 
= 3.178 pc cm~^, under the usual assumption of a 10% ISM ionization (where N^ is the 
electron column density). In other words, the PL model requires a very low ISM ionization, 
~0.4%, (i.e., a very large ratio Nn/N^ ~ 230, much larger than for any other radio pulsar 
detected in X-rays). Also, taking into account that B1929 is a nearby pulsar {d = 361 
pc), the A^H obtained from the PL fit is uncomfortably close to the HI c olumn density, 
A^HT = 3.5 X 10^^, throughout the entire Galaxy in the direction of B1929 (JKalberla et al. 



<LHL 



20051 ) . On the contrary, A^h 



0.17+o;^7 X 10^^ cm"2 inferred from the PL+BB fit is consistent 
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with the usual estimate based on the dispersion measure as well as with the iVn estimated 
from the PWN spectrum (see Fig. [TTj) . and it is much lower than the total Galactic A^hi- 

Also, the PL+BB model looks more attractive than the PL model because it gives the 
slope of the PL component, F g y 1.7, within th e range of spectral slopes found for a large sam- 
ple of young pulsars (F Ri 1 —2; lGotthelill2003l ). while the PL fit to the B1929 spectrum results 
in a considerably steeper slope, F Ri 3. One could argue that the PL fits to the spectra of 
other old pulsars (e.g. , B0950+08, B2224+65, B0823+26, B0628-28, B1133+16, B0943+10 



Zavlin fc Pavlov 



20061: 



2004 



Becker et al.ll2004J : iTepedelenhoglu fc Ogelman]l2005l : iKargaltsev et al. 



Zhang et al.l l2005l ) also show rather steep slopes, F fa 2 — 3, sugge sting that pulsar 



spect ra might soften with increasing age or decreasing spin-down power (IKargaltsev et al. 



20061 ). However, a plausible alternative interpretation of the softer spectra of old pulsars is 
that they are, in fact, comprised of a (soft) thermal compone nt and a PL component with a 
more gradual slopqj, similar to those found in young pulsars ( Zavlin fc Pavlovl 120041 ). 



Thus, based on the goodness of fit and the astrophysical arguments, we conclude that 
the X-ray emission from B1929 most likely includes both thermal and magnetospheric 
components. The BB model for the thermal component gives the apparent temperature 
kTa ~ 0.3 keV and projected emitting area A±^a ~ 3000 m^, at the distance of 361 pc, 
suggesting that the thermal emission originates from small heated spots (e.g., polar caps). 
These temperature and area correspond to the apparent radius R^ = (Ax.a/Ti')^^^ ~ 33 m 
and bolometric luminosity L' 



bol,a 



4^±,aO"^a ~ 1.1 X 10^° ergs s~^ of an equivalent sphere. 
The true size and luminosity of the polar caps depend on the geometry and the gravitational 
redshift factor, gr = (1 — Rg/RNsY^^, where Rg = 2.95Mns/M0 km, i?NS and Mns are 
the mass and radius of the neutron star. If there are two identical hot spots at the poles 
of a centered magnetic dipole, the polar cap radius and the luminosity of two polar caps, 
as measured at the neutron star surface, are Rpc = R'^f~^^^ and -Lboi,pc = L'^oia/C^fdr)^ 
where / < 1 is a geometrical factor depending on the angles ( beween the lin e of sight and 



the s pin axis and a between the spin and magnetic axes, as well as on gr (jPavlov et al. 
20071 ). For ins tance, assuming Mnr = l AMp^, i?NS = 10 km (i.e. gr = 0.766) , a = 36.0 
and C = 61. 5° (lEverett &: Weisberg|l200ll ) and using the approach described by IZavlin et al. 



(119951) and iBeloborodovl (J2002|), we find / = 0.897, which gives Rpc = IMR"^ ?a 35 m and 



^bol,pc 



1-62L- ^ 



1.8 X lO'^^ergss-i = 4.7 x 10"*^. 



The estimated polar cap radius is a factor of ~9 smaller t han the conve ntional polar 
cap radius of a radio pulsar, Rpc = (27ri?^g/cP)^/^ (see e.g., iMichell Il99ll ). which is ~ 



^Note that the observed spectra of these pulsars can be satisfactorily fitted with quite different models 
because these objects are very faint in the X-ray range. 
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300 m for B1929, assuming i?Ns = 10 km. However, the observed polar cap radii can be 
substantially smaller or larger than the conventioi ial value. In particular, BB fits of several 



old pulsar spectra show Rpc ~ (0.1-0.2) Rpc (e.g., IZavlin fc Pavlovl |2004J : IZhang et al.l 12005 
Kargaltsev et al.ll2006l ). For some pulsars, such a discrepancy can be alleviated assuming 
that the polar cap is covered by a hydrogen or helium atmosphere, in which case the effective 
temperature would be a factor of 2 lower, and the rad ius a factor of 3-10 larger while the 
bolometric luminosity would not change significantly (jPavlov et al.lll995l: IZavlin &: Pavlov 
20041 ). Another explanation for such a discrepancy was proposed by IZhang et al.l (120051 ) . 
who suggested that only a small fraction of the polar cap surface, associated with footprints 
of sparks produced by intermittent breakdowns of an "inner gap" above the polar cap, is 
hot enough to emit X-rays. 

The luminosity of the detected thermal component, which is less dependent on the as- 
sumed spectral model than t he radius and temperature, c an be compared with the models 
for pulsar polar cap heating (JHarding fc Muslimovl 120021 . l200ll ). For B1929, these models 



predict the polar cap luminosity Lboi.pc ~ 10'^^ ergss"^ if the polar cap is heated by positrons 
produced through curvature radiation of electrons accelerated in the neutron star magne- 
tosphere. Although the polar cap thermal luminosity est imated from the PL+BB sp ectral 
model is a factor of a few lower than that predicted by iHarding fc Muslimovl (120011 ). this 
can be considered as a reasonably good agreement, with account of the uncertainties of both 
the theoretical model and the observational results. We expect that deeper observations of 
other old pulsars would also detect thermal components from their polar caps and help us 
better understand the polar cap heating mechanisms. 

As Figure m shows, the thermal component dominates at 0.5 keV ^ E < 2keV. This 
means that we can expect different pulse shapes within and outside this energy range, because 
the pulsations of thermal radiation should be smoother and shallower than those of the 
magnetospheric radiation. The study of the energy dependence of pulsations not only can 
confirm the presence of the thermal component, but it would also allow one to infer the polar 
cap geometry and emission mechanism (e.g., BB vs. a l ight-element atmo s phere) , and even 
measure the mass-to-radius ratio for the neutron star (IZavlin et al.lll995l : iPavlov fc Zavlin 
19971 : IZavlin fc Pavlovl |2004| ). Since the current XMM-Newton data with sufficient time 
resolution do not have enough counts for such an analysis, a new deep observation of B1929 
is needed to perform this important test. 

The spectral slope of the PL (magnetospheric) component inferred from the PL-I-BB 
fit is r Ri 1.7, similar to those of younger pulsars. The luminosity of this component, 
Lx = 4:71 d'^F^'^^'^ ^ 1.7 X lO^'' ergs s~^ in the 0.3-10 keV band, correponds to the nonthermal 
X-ray efficiency r^nonth = Lx/ E ^ 4.4 x 10""^, typical for the whole population of radio pulsars 
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detected in the X-rays, and comparable to the thermal (polar cap) efficiency. Notice that 
the one-component PL fit results not only in a much softer spectrum (F ^ 3.0) but also in 
a higher luminosity, Lx ~ 6.6 x 10^° ergs s~\ and efficiency, r^nonth ~ 1-7 x 10~^, above the 
typical value for X-ray emitting radio pulsars. A similar trend has been observed in other old 
pulsars: PL fits of their spectra yield X-ray efficencies noticeably higher that those of young 
and middle-aged pulsars, which might suggest that the X-ray efficiency in this energy range 
grows with increasing age and decreasing spin-down power, perhaps at the expence of the 
efficiency at higher photon energies. However, this conclusion is based on the assumption 
that the contribution of the thermal component is negligible in the soft X-ray range, which 
we believe is not valid at least for the best-studied old pulsars, B1929 and B0950+08. To 
understand the evolution of X-ray properties of pulsars, deep observations of a larger sample 
of old pulsars are warranted. 



3.3. Summary 

Two ACIS-S observations of B1929 revealed a faint PWN surrounding this old, nearby 
pulsar. The observed morphology includes a compact nebula with two patchy wings in the 
immediate vicinity of B1929, and a tail extending in the direction opposite to the pulsar's 
proper motion, aligned with the much longer tail detected in the previous ROSAT and XMM- 
Newton observations. The shape of the nebula and its spectral properties are consistent with 
the proposed bow-shock classification. The properties of the compact nebula suggest that 
the pulsar wind outfiow is anisotropic, possibly concentrated toward the equatorial plane 
perpendicular to the pulsar velocity. The size of the observed tail implies an average fiow 
velocity of ~ 0.1c. The whole observed PWN radiates about 2 x 10^'^ of the pulsar spin-down 
power in X-rays. By comparing the observed PWN properties with the bow-shock MHD 
models, we estimated the temperature of the local ISM to be ~ 10^ K, which is consistent 
with the lack of the Ha bow-shock emission around this pulsar. 

In the spectral analysis of the combined Chandra and XMM-Newton data, we detected 
a thermal component in the pulsar radiation, whose luminosity, ~ (1^2) x 10^° ergs s^^, is 
comparable with the magnetospheric luminosity in the X-ray band. The thermal radiation 
is likely emanating from polar caps heated by positrons created by the curvature radiation of 
ultrarelativistic electrons accelerated in the pulsar magnet osphere. The spectrum and X-ray 
efficiency of the magnetospheric component are similar to those found in young and middle- 
aged pulsars. Further X-ray observations of the pulsar will allow one to better characterize 
its spectral and timing properties and understand the nature of the magnetospheric radiation 
and the mechanisms of polar cap heating. 
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Table 1: Chandra and XMM-Newton observations of B1929 



Instrument 


ObsID 


Date 


Exposure 


EPIC 


0113051301 


10 Nov 2003 


10.5 / 10.5 / 7.3 


EPIC 


0113051401 


27 Apr 2004 


15.2 / 16.4 / 11.0 


EPIC 


0113051501 


29 Apr 2004 


7.3 / 7.6 / 10.3 


ACIS 


6657 


4 Dec 2005 


20.9 


ACIS 


7230 


28 May 2006 


24.6 



Note. — The good-time exposures are given in ks. The first, second, and third EPIC exposures are for 
MOSl, M0S2, and PN, respectively. The effective PN exposure in each observation is ~70% of the total 
PN good-time exposure because of the reduced efficiency of the small window mode. 



o 



Table 2: Background-subtracted counts and surface brightness of PWN regions in Chandra 
observations 6657, 7230, and combined data 



Region 




1 (Southern Wing) 


2 (Northern Wing) 


3 (Front) 


4 (Inner Blob) 


5 (Outer Blob) 


Extraction area 
arcsec-^ 




31.5 


46.7 


77.9 


157 


1884 


Counts 
(bkg-subtracted) 


6657 
7230 
comb. 


18.9±4.6 
15.3±4.1 
34.1±6.1 


6.9±3.2 
17.3±4.5 
24.1±5.5 


12.7±4.2 
3.7±2.8 
16.4±5.0 


23.5±5.8 
14.4±4.8 
37.9±7.7 


57±14 
116±15 
173±20 


Bkg-subtracted 

surface briglitness 

(10~^ ctss"-*^ arcsec"^) 


6657 
7230 
comb. 


2.9±0.6 
2.0±0.5 
2.4±0.4 


0.7±0.3 
1.5±0.4 
1.4±0.2 


0.8±0.2 
0.2±0.1 
0.5±0.1 


0.7±0.2 
0.4±0.1 
0.6±0.1 


0.14±0.03 
0.25±0.03 
0.20±0.02 



hl^ 



Note. — The regions are shown in Fig. [3l The errors represent statistical uncertainties at the 68% 
confidence leveL 



Table 3: Fits to the spectrum of B1929 for the Chandra and XMM-Newton observations. 



Model 



AfH 



PL Norm. 



10- 



kT 

kcV 



Radius'' 



Xl/d.o.i 



Absorbed Flux 



Luminosity 
1030 



ergss 



Chandra: 

PL 

PL+BB 

XMM-Newton: 

PL 

PL+BB 

combined: 

PL 

PL+BB 



9 A^+0-52 
n 49+0-93 


2 94+0-25 
^■^^-0.22 
-I 09+1.04 


7.99ll;I? 

i-9ii?:33 


90+0.05 
U-^9_0.06 


-iq 0+46.3 


0.91/71 
0.88/69 


1 37+0-14 

l-^'-0.27 
1 c:9+0.23 
-^■'^^-0.37 


6.26; 
4.66: 




2 99+0-20 

1 fi'^+0-86 
i.Dd_g 29 


8.96+1-32 

1 4-3+3.73 
^•^■^-0.42 


.,„+0.02 


34.911190 


0.97/190 
0.94/188 


1 5:9+0.12 
i.oy_p 25 
1 09+0.13 


7.08: 
2.87: 


9 90+0.27 

n 17+0-23 
^■l'-0.17 


2.95+0-14 
1 70+0.46 


8-43ir79 
l-64t(V28 


n •Hn+0-02 

^■■''J-0.03 


33.1^^6 


1.03/264 
0.98/262 


1 51+0.12 
l-Sl-0.24 
-T 7t;+0.11 
l-'°-0.22 


6.62: 

2.84: 



+2.66 
-0.93 
+0.86 
-1.12 
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Note. — The observed flux and unabsorbcd luminosity, Lx = 4tt(PF^'^^^^, are calculated for the 0.3-8 keV 
energy band for Chandra, and 0.3-10 keV for XMM-Newton and combined data. The listed uncertainties 
are at a 90% confidence level determined for 2 interesting parameters. 
"Radius of equivalent sphere for the BB component (see $3.2). 



